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Abstract 

A possible interpretation of the near-threshold enhancement in the (pp)-mass spectrum in 
J/ip^jpp is the of existence of a narrow baryonium resonance X(1860). Mesonic decays of the 
(pp)-bound state X(1860) due to the nucleon-antinucleon annihilation are investigated in this 
paper. Mesonic coherent states with fixed G-parity and P-parity have been constructed . The 
Amado-Cannata-Dedoder-Locher-Shao formulation(Phys Rev Lett. 72, 970 (1994)) is extended 
to the decays of the X(1860). By this method, the branch-fraction ratios of Br(X — » n^ir), 
Br(X — > n2ir) and Br(X — » 3n) are calculated. It is shown that if the X(1860) is a bound state 
of (pp), the decay channel (X — > ry47r) is favored over {X — » n2ir). In this way, we develop crite- 
ria for distinguishing the baryonium interpretation for the near-threshold enhancement effects in 
(pp)-mass spectrum in J/ip^jpp from other possibilities. Experimental checks are expected. An 
intuitive picture for our results is discussed. 



* e-mail address: dinggj@mail.ustc.edu.cn 

t e-mail address: mlyan@ustc.edu.cn; corresponding author. 



1 



I. INTRODUCTION 



There is growing interest in exotic hadrons, which may open new windows for understand- 
ing the hadronic structures and QCD at low energy. Recently, the BES Collaboration 
observed a near-threshold enhancement in the proton-antiproton (pp) mass spectrum from 
the radiative decay J /if) — > ^pp This enhancement can be fitted with either an S- or 
P-wave Breit-Wigner resonance function. In the case of S'-wave fit, the peak mass is at 
M = 1859±? (stat)t^ 5 (sys) with a total width T < 30 MeV/c 2 at 90% percent confidence 
level. For the S— wave fit, the corresponding spin and parity are J PC = l ". Moreover, 
the Belle Collaboration also reported similar observations of the decays B + — > K + pp Q] 
and B° — > D°pp |3j], showing enhancements in the pp invariant mass distributions near 2m p . 
These observations could be naively interpreted as signals for baryonium pp bound states 
El The BES-datum fit in ref.[| represents the simplest interpretation of the experi- 
mental results as an indication of a baryonium resonance. Here, we denote this baryonium 
particle as X(1860) with J PC = h and I G = + [3]. However, this is only one possible 
interpretation. Otherpossible ways to understand this phenomena include, for instance, a 
flavorless gluon state js], a final state interactions 0] or an effect of the quark fragmentation 
process 0| , etc. In order to ascertain whether or not the X(1860) exists, more evidence is 
needed. A significant distinguishing feature for the baryonium interpretation is that the 
decays of X(1860) are mainly due to proton-antiproton annihilation in the baryonium. In 
this paper, we investigate pp-annihilations in Xf I860) by means of a coherent-state method 
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on nucleon-antinucleon annihilation in large N c QCD 10] 

At first glance, the most favorable X(1860) decay channel would be X(1860) — > r]7T7r, 
because it is the simplest hadronic process with the largest phase space. However, since 
the decay is caused by (pp) -annihilations in the X(1860), this naive observation may be 
not true. Low-energy nucleon-antinucleon annihilation is a fertile area for studying low 
energy QCD, and there are many experimental and theoretical studies in the literaturefe.g., 

see QHQQQ QQQ)- 

In these reports, it has been shown that the nucleon- 
antinucleon annihilation at rest mostly favors processes with between 4 and 7 pion final 
states, over those with two or three pion's Q|. This is a a general characteristic of (pp)- 
annihilation (without consideration of J PC and I G quantum numbers). It is interesting to 
pursue whether or not there are similar features in X(1860) decay. If so, we will have a new 
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criteria to characterize the X(1860). This is the main aim of the work reported here. 

The contents of this paper are organized as follows: In section II, we use a toy model 
to describe the possibility of the {p — p)-collisions inside a (pp)-bound state; in section III, 
we construct the mesonic coherent states with fixed G- and P-parities; Section IV describes 
calculation of the branch fractions of Br(X — > rjAir), Br(X — > i]2tt) and Br(X — > 3rj). 
Finally, we briefly summarize our results, and provide an intuitive picture for our results. 

II. A TOY MODEL DESCRIPTION ON (p - p) COLLISIONS INSIDE A (pp)- 
BOUND STATE 

In order to understand the mechanics of X(1860) decays, we assume that the proton and 
the antiproton collide with each other inside the X(1860), resulting in the collapse of the 
(pp)-bound state , or the decay due to the effect of rapid nucleon-antinucleon annihilation. 
We use the (p — p)-collision frequency v (i.e., the collision times per time-unit period) in the 
X(1860) to characterize the possibility of such (pp)-collisions. We show that this frequency 
can be estimated in a simple toy model, in which the proton and the antiproton are treated as 
point-like particles. This collision frequency is actually equivalent to the total X(1860) decay 
width in the annihilation- assumption mentioned above. Following ref. we roughly sketch 
the toy model for the X(1860) (pp)-bound state. The single-well potential toy model for the 
(pn)-molecular bound state deuteron first appeared in the literature more than fifty years 
ago [18]. The quasi-stable (pp)-molecular bound state X(1860) can be similarly described 
by a double-well potentiator Skyrmion-type potential) toy model[^|0. The potential of 
such a double- well model, V(p), is expressed as follows p] (see Fig.l) 

V(p) = 2m p -c8(p)+V c (p), (1) 

where 

h = m p /4:, < p < A, 

V c (p) = | -Vpp = -73 MeV, A < p < a pP , (2) 
0, p> a pp . 

where a p p ~ 2 fm, c ~ 2.161 , m p is the mass of proton and A = l/(2m p ) ~ O.lfm. In this 
case, the Schrodinger equation for S-wave bound states is 

-1 d 2 



2{m p /2) dp 



u(p) + [V(p)-E]u( P ) = 0, (3) 



where u(p) = p ip{p) is the radial wave function, and m p /2 is the nucleon reduced mass. 
This equation can be solved analytically, and has a bound state u^ip) with binding energy 
E 2 — —17.2 MeV due to the attractive square well potential at intermediate ranges (see 
Appendix). This molecular state is identified as the X(1860). 
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FIG. 1: The double well potential of pj?-system. 

We use this solution to estimate the (pp) -collision frequency inside the X(1860). Because of 
the effect of (pp)-annihilation, this bound state is not stable. Note that in this model, there 
are two attractive potential wells: one is at p ~ and the other is at intermediate values; 
between them is a potential barrier. When p ~ 0, the proton and antiproton are in collision. 
Therefore, we can derive the (pp)-collision frequency in the inside of X(1860) by calculating 
the quantum tunnelling effect for u 2 (p) passing through the potential barrier. In the WKB- 
approximation, the tunnelling coefficient (i.e., barrier penetrability) is (see Schiff's book in 
ref.Q) 

r x 
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the time-period 9 for the particle's round trip is 
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v ' y V pp + E 2 

Thus, the state X(1860)'s (i.e.,u 2 (r)) lifetime is r = 6>T -1 , and v, the the so called pp- 
collision frequency inside the X(1860), is equal to T x , the total width of X(1860): 

1 
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2\Jm p (h - E 2 ) 



(6) 



Substituting E 2 



17.2 MeV, a pp = 2.0 fm into Eq.©, we obtain the prediction 



x 



15.5 MeV, 
4 



(7) 



which is compatible with the experimental data 

The corresponding (pp)-collision frequency inside the X(1860) is about 

v = Tx/h ^ 2.35 x 10 22 Hz. (8) 

Note that because the binding energy E2 is rather small (compared to 2m p ), the annihilation 
processes that cause X(1860) to be unstable occur nearly at rest. 



III. COHERENT STATES WITH FIXED G— AND P-PARITIES 



When the proton and antiproton collide, they will rapidly annihilate into mesons. The 



coherent state method in and |12| investigates (pp)-annihilation at rest without consid- 
eration of the P— and G— parity. In this paper we are concerned with annihilation inside the 
X(1860) where the (pp) state has I G (J PC ) = + (0 _+ ). For this state, the allowed processes 
are (X(1860) — > 77271; tjAtt; ■ ■ ■ i]2mr; 3i]) (we ignore processes that involve K-mesons), 
where the tt and rj are both pseudoscalars and the G-parities are negative and positive re- 
spectively. In this case, we have to introduce P— and G— parities into the previous coherent 
state description of mesons radiated in fppVannihilations jiol. I12I] . 

We construct coherent states with fixed four momentum and fixed isospin, and, also, 
with fixed G parity and P parity. Following the method in 12| . we first construct the 
field operator F that creates a it or rj at the space time point x and directed in the isospin 
direction n. 

F(x, h) = J rf 3 ke-^V(k)a{ c • n + J d 3 qe~^ (9) 



where k ■ x = k$t — k ■ x with k = yk + , q ■ x = q^t — q ■ x with q = y q 2 + m 2 , is 
the isospin-triplet creation operator, and 6^ is the isospin-singlet creation operator. From 
the G- an P-parities of the tt and rj, we have 

Ghr(x,t)& = -7r(x,t) 
G7j(x,t)Gt = r?(x,t) 

PTT(^t)P ] = -7T(-X,t) 

Pr/(x,t)pt = -7j(-x,t), (10) 
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where G, P are the unitary operators as follows 

P = exp[z| E( a k. J a -k J + b k b -k + a k/k J + 6 k 6 k)] 
j,k 

C = exp[z- E( a k,-l a k,i + a k,l a k,-i " a k,i a k,l ~ a k,-i a k,-l)] 



7T 



ex P[-^ E( a k,o a k,l + a k,o a k,-l " a k,i a k,o - 4,-i a k,o)]- ( n ) 
It is straight-forward to check the following equations 

GtfqGi = -fol q (12) 

where % = 1, 0, —1 corresponding to tt + , tt°, it~ . Under G transformation, F(x, n) becomes 

G(x, n) = GF(x, n)G f = - / d 3 kf(k)ak • n e~ lk ' x + J d 3 qe- iq - x g{q)b j q . (13) 

For simplicity we take /(— k) = /(k) and g(— q) = g(q), then the P transformation of 
F(x, n) is as follows 

PF(x, n)P f = - /" d 3 k/(k)ak. • n e - **' - J rf 3 qe~^ :r '^(q)6 q = -F{x , n), (14) 

with x' = (—x, t). Then the desired coherent state with fixed four-momentum, fixed isospin, 
and also with well-defined G parity(+) and P parity(-) is constructed as follows 

\K,I,I,) = I^^B e iK -*\f,g,x,h,2)Y* Iz (h), (15) 

where 

\f, g, x, n, 2) = [e F ^ + e G ^ - F(x, n) - G(x, n) - e^*'*** - e~ G{x (16) 
-F(x',n) -G(x',n)]|0). 

Here we have subtracted the states without a meson and with only one meson, since they 
violate the conservation of energy and momentum. The states defined in Eq. (jl5j) are or- 
thogonal 



(K,I,I Z \K',I',I' Z ) = 6\K-K')6 H ,6 IzI ,l(K), (17) 



where 1{K) is the normalization factor: 

1(K) j ^ 11 e^^JWtriP- ){4exp[ P/ (x)n • n + p g {x)\ + 4exp[- P/ (o;)n ■ n + p g {x)\ 

-Spg(x) - 4exp[-p / (x)n • n - p a {x)\ - 4 exp [p f (x) n • n - p g {x)\ - 8p g (x)}, (18) 



where 



p f (x) = f fPp\f(p)\ a e** 

Pg (x)= fd 3 q\g( q )\ 2 e^. (19) 



We use the expansion method developed in |TJj to calculate the normalization integral 

m+n>2; m is even, n is odd 



where 



and 



Tfl Yl 77~l 71 

I(K,m,n) = J 5\K - 2> - £ ?i) II ^P*I/(P*)I 2 II ^(q,)| 2 (21) 

8=1 j'=l i=l j'=l 

F(m,I) = J ^ * Y; iz (h)Y IIz (h)(h-hT l 

I > m ot I — m is odd 

I < m and / — m is even. 



(22) 



(m-7)!!(/+m+l)!! 

Note that the effect of phase space for the decay has been taken into account via the 
8 4 (K - YT=\ Pi - £"=i Qj) UiLi d 3 p { \ n™=i d 3 qj in function I(K, m, n). Each individual term 
with (m,n) in the sum Eq.(|2()jl represents the contribution from the decay channel whose 
final particles are mir plus nrj. For a fixed total energy, the sum must terminate. The 
coherent state naturally gives the result that only decays to even numbers of 7r's and odd 
numbers of 77's. The decays conserve P— parity and G— parity. The mean numbers of tt of 
isospin type i and 77 are given by 

= {K,I,I Z \ J d 3 k a+ .a k .\K,I,I 2 ) 

= J ^^^^y^(n)y 7A (n') (4p / (^)A;%{exp[p / (x)n'*.n + p fl (x)] 

- exp[-p / (x)n* ■ n + p g {x))} + 4p/(x)nJ*n<{exp[-p/(a;)n * ■ n - p g (x)) 
- exp[p f (x)n* -h- p g (x)]}) (23) 



and 



N " = l(K) {KJJzl I d3k 6 k 6 kl^ J ^> 

= T(K) I (2^4 ^4^° e^Y* h (n)Y IJz (^) (4p 9 (x){exp[ P/ (*)n * -n + Pg (x)] 
+ exp[-pf(x)h * ■ n + p g {x)]} + 4:p g (x){exp[-pf(x)h'* ■ n - p a {x)\ 
+ exp[p f (x)n* -n- p g (x)]}) . (24) 



Using the expansion method, we obtain 
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N ** = rtr\ E - r AK,m+ 1^)0^,1,1.) (25) 
2{K) ^ mini 

m is odd 



and 



n is odd 



N v = ^ E -^-J(K,m,n + l)F(m,I), (26) 
T(K) ^ mini 

m is even 
n is even 



where 



\ n Y^Y^ih )n/n,(n • n*) m 
2/ -I- 1 

= E^K02jxi(<l' . 10 l /0 ><' n . 1 »1^>) 2 - ( 27 ) 

Zn 

IV. X(1860)-DECAY THROUGH (p - ^-ANNIHILATION 

Now we illustrate the X(1860)-decays due to (p — ^-annihilations. In Section II, we have 
shown how the p and p meet together in the X(1860) by using a model that can be use 
to interpret the near-threshold enhancement in the p — p mass spectrum in J ftp — > 7pp[lj]. 
The p — p-collision (or overlap ) must lead to the (p — ^-annihilation, and this causes 
X(1860)-decay. Moreover, in the previous section, the meson coherent states with fixed Gr- 
and P-parities are constructed that describe the final meson states radiated by the p — p- 
annihilation. After these preliminaries, the investigation of X(1860) decay becomes feasible. 

For the case of X(1860), 1 = 0, and from Eq.(|27jl we have 

G 1 (m, 0, 0) = G {m, 0, 0) = GU(m, 0, 0) = F(m, l)/3. (28) 



Instituting Eq.flU into Eq.(j23J, we get 

N w+ : N n o : N w - =1:1:1. (29) 

This indicates that among the products of X(1860)-meson decays, the ratios between the 
number of 7r + and ir~ , and between that for 7T ± and 7r° are fixed. 

The probability of the decay with annihilation-products of (N^n, N^rj) is given by 

P(N W , N v ) = — -— I Ud^Yld 3 ^ Kp 1 p 2 ---p J v ir q 1 q 2 ---q 2 vJ^,/,/z>r 

1" V i=l j=l 

1 16I(K,N n ,N v )F(N„,I) 
' X(7Q iVJiV,! • 1 ) 

In the same spirit, for the case of (N 7T +tt + , N n - n~ , N^tt , N v rj) annihilation-products , the 
probability is 

p(iv. + ,AUA^,Ag = — iv^IiMiv^ ' (31) 

where 

^ n r; /z (n)r Jj2 (n )(n + ■ n + *)^+(n • n *)^°(n_ • n!)^-. 

(32) 

Since the branching fraction 7?r (X(1860) — * run + 77,77) is proportional to P(m,n), from 
Eq.(jnnj) we can obtain the ratio of Br$(X — > ??47r) and Btq(X — > ^27r) as 

£?r (X -> 7747T) = 7(AT,4,1)F(4,0) 2! = 7(7^, 4,1) 

Br (X^7]27r) 4! 7(7T, 2, 1)F(2, 0) 20 7 (AT, 2, 1) 

and 



(33) 



7jr (X -> t^tt) = 7(K,2,1)F(2,0) 3! = 7(^,2,1) 

£r (X^3r/) 2! 7(7sT, 0, 3)F(0, 0) 7(7T,0,3)' 1 1 

Note that in the above coherent state calculation of 7?ro(X(1860) — > m7r + n.77), charge 
conservation was not taken into account. Here we expand our discussion to include charge 
conservation. Using 7?r(X(1860) — > mir + n.77) to denote the corresponding branch fraction, 
we have 

Br(X -> r/2vr) = 5r(X -> r/Tr+vr^) + 5r(X -> r]2n°) 

Br(X -> r/47r) = Br(X -> t^tt^tT) + 5r(X -> ^tt+tt^tt ) + 5r(X -> r)4ir°). (35) 
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Consequently, the ratios between these branch fractions are 

Br(X -> r]2n) I(K, 2, 1) 2 
Br(X -> 3tt) ~ /(AT, 0,3) 3 
^(X^^tt) /(if, 4,1) 7 



(36) 



5r(X -> 772tt) /(AT, 2, 1) 300 

Since ~ 139MeV and ~ 547MeV, the X(1860) can decay to i]2ir, 3r], rjAir, i]67T, 
and T]8ir, but not into 3t]2tt and still conserve energy. We investigate the decay channels of 
X(1860) -»• (r/2vr, r/47r, and 3r?) below. 

Following ref.[ll|, we suppose that the meson field source turns on at t = and then 
decays exponentially in time, and that it has a spherical symmetric Yukawa shape. In this 
case, /(k) (as a Fourier transformation of the meson field source) is 

2 c k 2 

1/(^)1 = 7j 2 , oVtT 2 I TV) 2~' (^^) 

(k + a 2 ) 2 (uj^ + 'J 2 ) 2 uj^ 



where uj^ = y k 2 + m 2 , C is a strength and can be fixed by required that the average energy 
)e the energy released in annihilation, namely 2m p . In units of pion masses = 1), by 



11 1 we take a = 7 = 2. This corresponds to an annihilation region with a time and distance 
scale of half a pion Compton wave length — a reasonable size that gives reasonable agreement 
with experimental data 10, ll2j. Since both 7r and 77 belong to the pseudoscalar meson octet, 
we expect that g(k) should be the same as /(k) except which m n replaced by m^. With 
this parameter choice, without consider charge conservation, we can obtain from Eq. (j3*3*j) and 
Eq.flHU): 

Br (X -> 774tt) , 

— °-) '—L ~ 3.8 x 10 4 

Br (X -> 7]2ti) 

Br (X^3,) - 8 ' 8 - (38) 
When charge conservation is taken into account, using Eq.()36|). the results become as follows 

Br{ , X - ***\ M 1.8 x 10" (39) 

^ " 2 *> « 5.9. (40) 
5r(X -> 3r/) V ; 

From Eq.([39|) we find out that £>r(X — > ?727r) is heavily suppressed by about 4 orders 
of magnitude compared to BriX — > tjAtt). Namely, the most favorable decay channel is 
(X — > rjAir), rather than (X — > r]2rr). This prediction may not be quantitatively exact, but 

10 



must be qualitatively correct. We expect this to be a significant feature of the decays of (pp)- 
bound states with I G (J PC ) = + (0 h ) due to the nucleon-antinucleon annihilation decay 
mechanism, which is significantly different from the naive argument for the decays of ordinary 
particle as discussed above in the Introduction. Thus, we conclude that an experimental 
check of this prediction is meaningful for distinguishing the baryonium interpretation for the 
near-threshold enhancement effects in (pp)-mass spectrum in J/ip^'ypp from other possible 
interpretations. In the other hand, such experiments should be valuable efforts also because 
they will belong to seek new evidence for the existence of exotic hadron X(1860). 

Equation (|4()jl shows Br(X — > 3rj) << Br(X — > t]2tx). This is mainly due to the effects 
of the decay phase space, and the result is reasonable. 



V. SUMMARY AND DISCUSSION 



One of the possible interpretation of the near-threshold enhancement in the (pp)-mass 
spectrum in J/i/j^^pp is the existence of a narrow baryonium resonance X(1860). The 
mesonic decays of the X(1860) due to the nucleon-antinucleon annihilation have been inves- 
tigated in this paper. In order to clarify the picture of proton-antiproton annihilations inside 
X(1860), we employed a toy double well potential to derive the (pp) -collision frequency, (or 
collision possibility) inside the X(1860). In this model, the annihilations cause X(1860) de- 
cays. Specifically, in the model the proton and the antiproton are separated by a potential 
barrier, and the (p — p)-collision's frequency is computed by considering quantum tunnelling 
effects. We further construct meson coherent states with fixed G-parity and P-parity. These 
enable us to extend the Amado-Cannata-Dedoder-Locher-Shao formulation to discuss the 
decays of the (pp)-bound state X(1860). In this formalism, the process of pseudoscalar 
meson radiation from the annihilation is rapid, and, hence, the mesons are classical, and 
can be approximately described by coherent states. By this method, the ratios between the 
Br(X — > r]4ir), Br(X — > i]2tc) and Br(X —>■ 3r)) branching fractions are derived. Taking 
appropriate meson field source functions, and evaluating the integrals related to 3-body and 
5-body phase space in the decay processes, we obtain quantitative predictions. We find that 
in contrast to naive arguments, the Br(X — > r/27r) is heavily suppressed about four orders of 
magnitude in comparison to Br(X — > rjAn). In other words, if X(1860) is a bound state of 
(pp), the most favorable decay channel must be (X i]4ir), rather than (X — > i]2tt). This 
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provides a criteria for distinguishing the baryonium interpretation for the near-threshold en- 
hancement effects in (pp)-mass spectrum in J/tp—^'-fpp from other possibilities. Experimental 
checks are needed. 

The unexpected result in Eq. (J39|) results from calculations based the coherent state theory 
that successfully describes nucleon-antinucleon annihilations. This can be seen from an 
intuitive picture. Naively, the number of valence quarks in X(1860)(or (pp)) is equal to the 
number of valence quarks in (rjinr), i.e., in both systems there are three quarks plus three 
anti-quarks, so it seems that the decay X — > r]27i should be the most easily accomplished. 
However, the gluon content for (pp) and (r)7nr) are different. Generically, the gluon mass- 
percentage in the proton (or antiproton) is larger than that for the 7r or?^. This point can 
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20, 



2lj. In the 



be seen from some QCD-inspired models, for example the Skyrme model 
chiral limit of this model (i.e., m quarc k = 0), the masses of the baryons, including the proton, 
are non-zero, but the masses of 7r and i] 8 (that is the main component of if) vanish. Thus, 
one could interpret this as being due to more gluons in baryons which make them massive 
even in the limit of massless quarks. This indicates that there are some "redundant gluons" 
that are left over in the process of X = (pp) — > t]2tt. Consequently, the process might be 
expressed as 

X = (pp) -> 7/2ttG, (41) 

where G represents the "redundant gluons". Most likely, G and r] could combine to form 
the meson ?/, in which case the process (14 lj) becomes 
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X -> 2n(r]G) = 2nr}' 

I > 777T7T (42) 

Eg . (J42j) is just (X — > ?y47r), where the factor of 77' — > ?77r7r is the dominate channel to be 
considered (i.e., Brirj — > r]7m) ~ 65%)[24j. In this view, the process of (X — > ?]27r) will be 
almost forbidden and the (X — > tjAtt) or (X — > r]'2ir) would be most favorable, i.e., 

Br(X -> r]'2n) » Br(X -> t^tt). (43) 

Since >> m^, this is a very unusual result. This can be tested experimentally. 

This analysis could be extended to a description based on classical SU(3)f fields to 
describe the small X(1860) decay branching fractions into K mesons. 
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APPENDIX A 



We solve the Schrodinger equation (jHJ) with the toy model's potential ((TJ) (see Fig. 1) in 
the text. Namely, the potential V(p) ( note p is at r p = 0, p is at r p = xi + yj + zk, and 



p = \/x' 2 + y 2 + z 2 ) is 



V(p) = 2m p -c5(p) + V c {p), 



where 



V c (p) 



h = m p /A 



-73 MeV 



< p < A 
A < p < a p p 
p > a pP 



(Al) 



(A2) 
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where A = l/(2m p ) = 0.1 fm , a pp = 2.0 fm, and the equation is 

' d \u(p) + (V(p) - E)u(p) = 0, 



(A3) 



2{m p /2) dp 2 

where u(p) = pnp{p) is the radial wave function, m p /2 is the reduced mass, and < p < 
oo. Equation ()A3|) is a one dimensional wave equation with both one-dimensional delta- 
function potential (—c5(p)) and square- well potential, and can be solved analytically. We 
are interested in the bound state solutions. The corresponding wave function boundary 
condition is: 



u(p 



oo 



0. 



(A4) 



A mathematic trick for solving Eq. (|A3J) is as follows: We can mathematically extend the 
variable p from the region (0 < p < oo) to (— oo < p < oo) with setting V(— p) = V(p). In 
this way, —c5(p) in the V(p) becomes regular. We then can find out the solutions u(p) by 
standard procedure to solve Eq. (|A3|) with boundary condition u(\p\ — > oo) — > . Finally, 
we take u(p > 0) (noting, in the physics region p = y/x 2 + y 2 + z 2 > ) as the physical 
solutions which satisfy both differential equation eq. ()A3j) and the boundary condition (|A4|) . 

There are two bound states U\{p) and U2(p)- U\{p) with binding energy E x < —V p p = 
—73 MeV is due to — c5(p)-function potential mainly, and U2(p) with binding energy E% > 
—73 MeV is due to the attractive square well potential at middle range mainly. They are 
as follows 



^(g-v/mpMi) p -|- aie \/™pMi) P\ 
A 2 ( e -\/ m p(-Vpp- E i) p _|_ / 5 ie \/ m p(- 
A 3 e~V-^ p 



V pp — Ei) p\ 



Q 1 ^ e -y/m P {h-E 2 ) p _|_ a2e y/m p (h-E 2 



)P\ 



B 2 sm(^m p (Vpp + E 2 ) p + fa), 



< p < A 

A < p < a pp , (A5) 
P > a P p 

< p < A 

A < p < a pp , (A6) 
P > a pp 



where ax, a 2 , fix, fti, and At, Bi with i = 1,2,3 are constants. Due to the potential of 
— cS(p) in the V(p), there are relations between c and ax, a 2 as follows 



2(1 -a 2 ) 



h-E 2 2(1 -ax) 



h-E x 



«2 \ 



rrir. 



oix \ 



(A7) 
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We return to solving the quantum mechanics problem given above. For u\{p) or u 2 (p), 
the wave-function continuum conditions £1X6 £IS follows 



T P logUi{p) ^ = J P 

where i — 1, 2. Then we have 



— log Ui(p)\( p=x -) = —\ogu i (p)\ (p=x+) , (A8) 



—JogUi(p)\ {p=a -, = —\og Ul (p)\ {p=<) , (A9) 



m p (h - Ei] 



( e -\fmp{h~E{) A aieX /m p (ft-£;i) A^ 

f _g- \/ m P (-^pp-^i) A _|_ ^ ie A/^p(-Vpp-Si) A-j 



m P (-VL- - £i) ^ ^ . ^ , (A10) 

(e~V "ipC-Vpp-Si) A _|_ ^ e \Jm p {-V P p-E{) A^ 
^_ e -V^pHW-Ei) «pp _|_ p ie V mpi-Vpp-E^ Opp\ 



m p (-Vpp - Ei) 



^-^/mpi-Vpp-Ei) a PP _|_ p i6 \/ m p(-Vpp-Ei) a pP ^ 

-^-m p Ei , (All) 

(_ e ~\/mp(h-E2) A _|_ Q, 2e \/mp(/i-E2) A-j 



m p (/i - E 2 ) 



mpiVpp + E 2 ) cot Jm p (V pp + E 2 )\ + f3 2 ) , (A12) 



^ e -y / mp{h-E 2 ) A a2eV /mpCi-S 2 ) A\ 



m p (Vpp + £ 2 ) cot Jm p (Vpp + E 2 ) a pp + 



-J-m p E 2 . (A13) 



These equations can be solved numerically. E\ is an input of the model. Following ref. 
and taking E\ ~ — 976MeV << — 73MeV as input (see below), we get the parameters in the 
solutions as follows 

ai = 0.025, a 2 = -0.352, (A14) 
p 1 ~ 0, (3 2 = -3.383, 
c = 2.161, 

from which we get the binding energy of the state u 2 (p), 

E 2 = -17.2 MeV. (A15) 
This is the result that is used in the text. 
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The functions u\(p) and U2{p) are shown in Fig. 2 and Fig. 3 respectively. From the figures, 
one can see that u\(p) is a sharply peaked curve with a maximum at p = 0, and U2(p) has 
a node at p ~ 0.3fm and |w 2 (/o) | has an absolute maximum at p ~ 1.5fm. Both Ui(p) and 
u 2 (p) satisfy the boundary condition given by Eq. (jA4|) . 



cl 25 




FIG. 2: The Wave Function m(p). 
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p(fm) 

FIG. 3: The Wave Function U2{p)- 
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In order to be sure of that the three-dimensional wave functions of our solutions can be 
normalized, we verify that 

p poo 1 p roo 

rf 3 x|^(x)| 2 = J dQj o p 2 dp-\u i {p)\ 2 = J dQ J o dp\ Ul (p)\ 2 (A16) 

/ fctpp roo \ 

= 47T / d P \ Ul (p)\ 2 + dp\u t {p)\ 2 \ 



An (const. + const, f dpe~ 2 ^~ nipEl J 



= finite, 

where i = 1,2 and Ui(p) have been given in eqs. (EU), dS|. Equation ijXTSj) is a check to 
the rationality of our solutions. 

In the following, we discuss the the parameters in the model in order: 

1. a pp and V pp : We take the width of the square well potential, denoted 

close to that of the deuteron, i.e., a pp ~ a pn ~ 2.0 fm. According to QCD inspired 
considerations P, 0, Q the well potential between q and q should be twice as 
attractive as the q — g-case, i.e., the depth of the pp square well potential is V pp ~ 
2V pn = 73 MeV. 

2. h and A : The quantitative results of the model somehow depend on the parameters of 
the potential, such as the barrier height 2m p + h and width A. For definiteness, we have 
taken them to be about 2m p + m p /A and l/(2m p ) in solving the Schrodinger equation. 
Actually, this is reasonable because the dependence on height and width are weak in 
the practice calculation. The results with several different values for height and width 
are listed in Tables 1 and 2, where we see that the above values can give a reasonable 
binding energy E2 and decay width Fx, in compatible with the BES measurement. 



barrier height 2m p + h 


2m p + m p /8 


2m p + m p j 6 


2m p + m p / 4 


2m p + m p /2 


2m p + m p 


E2 


-17.2 MeV 


-17.2 MeV 


-17.2 MeV 


-17.3 MeV 


-17 A MeV 


T x 


17.7 MeV 


16.8 MeV 


15.5 MeV 


12.7 MeV 


9.7 MeV 



TABLE 1: The binding energy E 2 and width T obtained 
by solving the Skyrmion-type potential model with the 
potential barrier height 2m p + h from 2m p + m p /8 to 
2m p + m p . 
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barrier width A 


0.5/(2m p ) 


l/(2m p ) 


2/(2m p ) 


3/(2m p ) 


4/(2m p ) 




-17.1 MeV 


-17.2 MeV 


-17.4 MeV 


-17.4 MeV 


-16.9 MeV 


T x 


19.3 MeV 


15.5 MeV 


9.9 MeV 


6.4 MeV 


4.2 MeV 



TABLE 2: The binding energy E^ and width T obtained 
by solving the Skyrmion-type potential model with the 
potential barrier width A from 0.5/(2m p ) to 4/(2m p ) with 
h = m p /4: fixed. 



3. c and Ei : At p ~ 0, V(p) ~ — c 5(p) with a constant c > 0, which is a free 
parameter in the model. E\ is the eigenvalue of Ui(p), which is roughly the energy 
level of the Schrodinger equation with the one dimensional delta-function potential 
V(p) ~ — c 5(p). So, Ei is c-dependent and, hence, once E x were fixed, the free 
parameter c fixed. For definiteness, in this paper, we have taken 

E x = -(2m p - m m ) ~ -976 MeV, (A17) 

and then the corresponding c value is c ~ 2.161 (see Eq. (jA14jl ). 
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